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Abstract 

We prove that the “Immediate Exchange Model” of [5] has a dis¬ 
crete dual, where the duality functions are natural polynomials asso- 
caited to the Gamma distribution with shape parameter 2 and are 
exactly those connecting the BEP{2) and the SIP{2) models in [T], 

m- 

As a consequence, we recover invariance of products of Gamma distri¬ 
butions with shape parameter 2, and obtain ergodicity results. Next 
we show similar properties for a more general model, where the ex¬ 
change fraction is Beta{s,t) distributed, and product measures with 
Gamma(s-l-t) marginals are invariant. We also show that the discrete 
dual model is itself self-dual and has a similar continuous model as its 
scaling limit. We show that the self-duality is linked with an under¬ 
lying SU{1, 1) symmetry, similar to the one found before for the SIP 
and related processes. 


1 Introduction 

Kinetic wealth exchange models (KWEMs) constitute a popular class of 
econophysical models in which agents exchange their wealth according to 
some stochastic rules, always preserving the total amount of wealth in the 
economy. The aim is to understand some important properties of the dynam¬ 
ics of wealth distribution, such as wealth concentration, stationary distribu¬ 
tions and time dependent correlation functions. For a recent review about 
KWEMs, we refer to [3]. The apparently economically strong assumption 
of wealth conservation - which also rules out the possibility of (endogenous) 
growth - is justifiable by choosing the appropriate time scale (or time unit) 
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for the economy. An interesting feature of KWEMs is their similarity with 
another family of models, known as KMP [T]. Introduced in [7j, KMP mod¬ 
els are microscopic models of heat conduction and are meant to provide a 
microscopic foundation of the Fourier law; in those models the exchanged 
quantity represents energy. As shown in [T], duality is a powerful tool to 
study the properties of KMP models. Thanks to duality it is possible to 
investigate invariant measures, ergodic results, and important macroscopic 
properties such as hydrodynamic limits, the propagation of local equilib¬ 
rium, and the local equilibrium of boundary-driven non-equilibrium states. 
In [1] , the authors show that duality can also be fruitfully applied to kinetic 
wealth exchange models, obtaining relevant information about the station¬ 
ary distributions of a model with saving propensities. 

In this paper we aim to extend the use of duality techniques in the field 
of KWEMs, by focusing our attention on a recent model, the so-called “Im¬ 
mediate Exchange Model”. 

The model has been first proposed in [5], where it is studied via simula¬ 
tions, and it has been later analytically explored in |6]. In that model, 
upon exchange, each agent gives a fraction of his/her wealth to the other. 
In [6] it is proved that, if this fraction is a uniformly distributed random 
variable with support [0,1], then the exchange process is characterized by 
an invariant measure, which can be expressed as the product of Gamma{2) 
distributions. 

It is now worth noticing that an invariant measure which is a product of 
Gammas also occur in the redistribution models presented in [T]. In these 
models duality is characterized by duality polynomials that are naturally 
associated with the Gamma distribution and it is shown that these polyno¬ 
mials are also the duality functions linking a discrete particle system, the 
symmetric inclusion process SIP{k), with a diffusion process, the Brownian 
energy process BEP{k). It is therefore natural to conjecture that these poly¬ 
nomials also occur as duality functions in the Immediate Exchange Model 
of [5], relating this model to a simpler discrete dual model. In this pa¬ 
per we show that this is indeed the case, and we generalize the Immediate 
Exchange Model to the case in which the random fraction of wealth the 
agents exchange is Beta{s, t) distributed. In this more general setting, the 
invariant measure shows to be a product of Gamma{s -|- t) distributions. 
As in [3], using duality we are able to directly infer basic properties of the 
time-dependent expected wealth, together with an ergodic result. 

The rest of our paper is organized as follows: in Section 2 we describe 
the Immediate Exchange Model when the economy is just made up of two 
agents and prove duality with a discrete two-agent model. In Section 3 we 
extend the model to the case of many agents and we give some relevant con¬ 
sequences of duality. In Section 4 a further generalization is proposed, by 
assuming Beta{s, t)-distributed exchanged fractions of wealth; also for this 
generalized model we obtain duality with a discrete model and stationary 
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product measures which are Gamma with shape parameter s + t. In Section 
5 we study various properties of the discrete dual process, which is an inter¬ 
esting model in itself. We characterize its reversible product measures and 
prove that in an appropriate scaling limit it scales to a simple variation of 
the original continuum model. Finally, in section 6 we show self-duality of 
the discrete model for the general case via a Lie algebraic approach, where 
we actually obtain the full SU{1, 1) symmetry of the discrete model, and, as 
a further consequence, of the continuous model, too. Self-duality then fol¬ 
lows by acting with an appropriate symmetry on the so-called cheap duality 
function obtained from the reversible product measure [2]. 


2 The Immediate Exchange model with two agents 
and its dual 


2.1 Definition of the model 


We start by considering a toy economy with just two agents, as given in 
[5] and [6]. More complex models can be built by addition of two-agent 
generators along the edges of a graph. Most properties such as duality and 
self-duality transfer immediately from the two-agent model to the many 
agent models. 

More formally, we write {x,y) G Q, with 12 = [0,oo)^. With s = x + y 
we indicate the total wealth in the economy. 

Then the dynamics of two agents is described as follows, starting from an 
initial state {Xq,Yq) = {x,y), after an exponential waiting time (with mean 
one), an exchange of wealth occurs, whereby the wealth configuration (x,y) 
is updated towards {x',y'), with 

x' = x(l — U) + yV 

y' = y(l-F)+xI/, (1) 


where U and V are two i.i.d. Uniform{0, 1) random variables. This gives 
a continuous-time Markov jump process {Xt,Yt) for which the total wealth 
Xt + Yt = Xq + YQ = x + yis conserved. 

The infinitesimal generator of this exchange process is defined on bounded 
continuous functions / via 


Lf{x, y) = lirn 
i—>-0 


¥.,^yf{Xt,Yt)-f{x,y) 
t 


1 

(/(x(l -u) + yv, y{l - u) -1 xu) - f{x, y)) dudv. (2) 

Notice that L can be rewritten as P—I, where P is the discrete-time Markov 
transition operator 



Pf{x,y) 



u) + yv,y{l 


v) -|- xu)dudv, 
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Figure 1: The continuous model and its discrete dual 
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and I is the identity. 

We denote {Xq, Yq) = {x, y) to be the initial wealth configuration of the two 
agents, and (Wj Yt) indicates the wealth of the two agents at time t > 0. 

2.2 Duality for the two-agent model 

We will now first define a discrete wealth distribution model, i.e., where 
wealth can only be a nonnegative integer quantity. See the figure for the 
continuous model and its discrete dual. This model will be related to the 
original one via a duality relation. 


Hence, in the discrete model the couple (x, y) G Q is replaced by a couple 
{n,m) G N^, where N denotes the set of non-negative integers (including 
zero). 

On this couple we dehne a continuous-time Markov process with generator 

n m - - 

^f{n, m) = Y] -k + l,m-l + k)- f{n, m)). (3) 


In this process, when initiated at (n, m), for a given k, I with 0 < k < n,0 < 
I < m, the wealth configuration changes from (re, rre) to {n—k+l,m — l+k) at 
'^ate • We denote this discrete state space continuous-time Markov 

process by {Nt,Mt), with {Nq,Mq) = {n,m). 

We will then show that the processes {Xt,Yt) and are related via 

duality. To introduce this, we need some further notation. 

Define, for x G [0,oo),re G N, the polynomial 


d{n, x) 


= X 


m 

r(2 + re) 


(re 1)! 


(4) 


and 


D{n, m\ X, y) = d{n, x)d{m, y). 


(5) 
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The d{n, •) polynomials are naturally associated to the Gamma distribution 
vg with shape parameter 2 and scale parameter 9, i.e. 


1/61 (dx) 


^xe-^/^dx 

9^ 


by 


J d{n,x)i^g{dx) = 0” 


for all n G N. 

With a slight abuse of notation, we will denote by i'g{dxdy) the product 
measure with marginals i^g. 

We are now ready to state the first main result. 


Theorem 2.1. The processes {Xt,Yt) and are each others dual 

with duality function given by ([^. 

More precisely, for all {x,y) G [0, oo)^, (n, m) G N^, and for all t > 0, we 
have 

Ea,,yT)(n, m; W, Yt) = ^n,mD{Nt, Mp, x, y), (6) 

where '^x,y o,nd '^n,m the expectations in the path-space measures started 
from {Xq,Yq) = {x,y) and {Nq,Mo) = {n,m) respectively. 


Proof. 

To prove Q it is sufficient to show the same relation at the level of the 
generators. In other words, we have to show that 


LD{n, m; x, y) = .^D{n, m; x, y), 


(7) 


for all (x,y) G [0,oo)^ and (n,m} G N^, and where L works on (x,y), and 
on (n, m). 

We compute 


PD(n, m; x, y) = 


1 



1 


/o Jo (n + l)!(m + l)! 
1 


(x(l — u) + yn)"(y(l — n) + ux)™' dudv 


n m / \ / \ /*1 /•! 

n\ (m\ / 


(n + l)!(m + 1)! 
1 

(n + l)!(m + 1)! 
1 


EE 

/ c =0 /=0 
n m 

EE 

/ c =0 /=0 
n m 

EE 


X- -y-y- "X" / / (1-u)^-V(l-nr-'u'dudu 

kj \ I J Jo Jo 

f ^n-k+iym-i+k k\{m-l)\ l\{n-k)\ 


k \ I 


ni 


{k + m — I l)\ {n — k 1 + l)\ 
m'. k\{m-l)\ l\{n-k)\ 


(n + l)!(m + 1)! ^ ^ (n — k)\k\ (m — l)U\ (k + m — I + 1) \ (n — k + I + 1)\ 
k=01=0 




D{n — k + l,m — I + k; x,y). 


■l+k 
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Now we have 


ft Hi ^ - 

EE—— = ' 


fc=0 1=0 

Therefore, we indeed find that 

n m 


n + 1 m + 1 


( 8 ) 


LD{n,m;x,y) = EErE:E — {D(n — k + l,m — I + k]x,y) — D{n, m; x, y)) 


n + 1 m + 1 

k=01=0 

= ^D{n,m] x,y). 


□ 


(9) 


As a consequence of duality, and thanks to the relation between the du¬ 
ality functions and the measure we obtain relevant information about 
the invariant measures. Let us denote by the set of probability mea¬ 
sures on [0, oo)^ with finite moments of all order and which are such that 
their finite moments determine the probability measure uniquely. I.e., two 
measures in with identical moments are equal. We say that such a 
measure satisfies the “finite moments condition”. Similarly for a probability 
measure on [0, oo) we say that it satisfies the “finite moments condition” 
if it has finite moments of all order and which are such that these finite 
moments determine the probability measure uniquely. This is e.g. assured 
by the Carleman’s moment growth condition. We will focus from now on 
only probability measures in this set ^f. 

Theorem 2.2. A probability measure v G is invariant if and only if its 
D-transform 

is harmonic for the dual process, i.e., if and only if 

En,mk{Nt, Mt) = v{n,m). 

for all n, m G N. In particular the product measures ve{dxdy) are invariant 
for the process {Xt,Yt). 

Proof. To have invariance of n G for all {n,m) G N^, it is sufficient 
to have 

j E:,^yD{n,m;Xt,Yt)i'{dxdy) = j D{n,m] x,y)n{dxdy) = t){n,m). (10) 
Combining this with duality and Fubini’s theorem we obtain 


v[n,m] = 


E,.^yD{n, m; Xt, Yt)n{dxdy) 
En,mD{Nt,Mt;x,y)iy{dxdy) = En,mk{Nt, Mt). 
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As a result, we find that v is invariant if and only if 

To show the invariance of the I'g measures, just notice that 

and recall that in the process (Nt, Mt) the sum Nt + Mt is conserved. □ 


Another consequence of duality is the ergodicity of the process {Xt,Yt). 
i.e., starting from any initial condition {x,y) the process converges to a 
unique stationary distribution determined by the conserved sum x + y. In¬ 
deed, the dual process starting from (n, m) is an irreducible continuous-time 
Markov chain on the finite set {{k, 1) : k+l = n+m} and therefore converges 
to a unique stationary distribution denoted by i^n+m, given by 


n. n _ 

^n+m(ki 0 ~ av ’ 

where 

^n+m = (k + 1)(? + 1). 

k,l:k-l-l=n-hm 

For all (n, m) G we can therefore observe that 


( 11 ) 

( 12 ) 


lim E.t 


iDiNt,Mt;x,y)) 

D{k,l-x,y)un+m{k,l) 


.x,yD{n,m]Xt,yt) = lira 

E^OO 

= E 

k,l:k-\-l=n-\-m 

It then follows from an easy computation using © that 

(x + 7/)"+"* 


E 

k,l:k-\-l=n-\-m 


D{k,l;x,y)un+mik,l) = 


[n + m)\Xi+m' 


(13) 


(14) 


where 3fn+m is given by (12) i.e., the limit in the r.h.s. of (13) only depends 
on x+y. On the other hand, in the process {Xt, Yt) we know that S = x+y is 
conserved. Therefore, the conditional measure obtained by conditioning the 
stationary product measure vg on the sum being equal to s is an invariant 
measure concentrating on the set {{u,v) G [0,oo)^ : u -|- u = s}. This 
measure is exactly the distribution of (se, s(l — e)), with e being Beta{2, 2). 
If we combine this fact with (13), we obtain the following ergodic theorem 
and complete characterization of the set of invariant measures satisfying the 
finite moment condition. 


Theorem 2.3. a) The process {Xt,Yt) is ergodic, i.e., {Xt,Yt) converges 
in distribution to {Se, S{1 — e)), with e ^ Beta{2,2) andS = XQ + Yo. 


b) The set of invariant measures contained in is given by the distri¬ 
butions of the form {Se, 5(1 — e)) where S is an arbitrary distribution 
on [ 0 , 00 ) satisfying the finite moments condition and e~ Beta{2,2). 
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3 Generalization to many agents 


Consider now an economy populated by many agents. Let us assume that 
the economy can be represented as a graph S, where each vertex represents 
an agent. Consider now an irreducible symmetric random walk kernel p{i,j) 
on S, i.e., such that p{i,j) = p{j,i) > pih j) = 1, and for all i,jGS 
there exists n with p^^\i,j) > 0. 

In this setting, the wealth configuration of the economy is an element of the 
set n = [0,oo)‘^. For x G O (from now on simply x), we denote with Xi the 
wealth of the agent i, that is of vertex i. 

We then define the generator of the model via 

Lf{x,y) = Y^p{i,j)Lijf{x), ( 15 ) 

ij 

with 

Lijfix) = j (/(x*^’“’') - f{x)) dudu, 

where 

' Xk 

^ y..y jf ^ ^ 

^Xj{l — u) + XiU li k = j 

Accordingly, the dual process has state space and the elements of this 
state space are denoted by ^ (from now on just ^), where is the number 
of “dual units” at vertex i. A configuration ^ is called finite if |^| = is 

finite. 

The generator of the dual process is then 

= Y^p{i,j)Lijf{n), ( 16 ) 

ij 

with 

r«/(0 = E E + + 

K=0L=0^^^ ASj I J 

where 

iik 

Ci'^^ =hi-K + L \ik = i 
^ij-L + K \ik = j 

Now, for ^ G N'^ and x G fl, define 

D{^,x) = ]J(i(^i,Xi) 

ies 
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(17) 



The relation between these duality polynomials and the product measure 
ue := ®i(zsi>e{<^Xi) is 

J D{^,x)ug{dx) = 9^^^ (18) 

with 

ifi = E«* 

ies 

the number of dual particles. 

In the many agents economy model, the duality relation between both 
processes is then given by the following theorem. Its proof is direct from the 
two agents case, because the generator is a sum of two agents generators. 

Theorem 3.1. Let ^ G N'^ be a finite configuration. For all x £ Ll and for 
all t > 0, we have 

Ea;T)(^, xt) = x). (19) 

As a consequence, the product measures vg = ®i^si^e{dxi) are invariant. 


Notice that when S is finite, the product measures ( 8 )jg 5 r' 0 (dxi) can never 
be ergodic because the total wealth is conserved. However, for infinite S, 
we have ergodicity under an additional condition. Let us denote by f) 
the probability to go from the finite configuration ^ G to the finite 
configuration in time t > 0, in the dual process with generator (15). 
Assume that 


lim pt(^,^') = 0 


( 20 ) 


for all G As an example we have S' = Z'^ and p{i,j) symmetric 
nearest neighbor random walk. 


Proposition 3.1. Let S be infinite and letp{i,j) be such that (20) holds. 
Then the product measure ®i^si^e{dxi) is ergodic 


Proof. Abbreviate v := ®i^s^e{dxi). Because ergodicity is implied by 
mixing, it suffices to show that 


^lim J KxD{f^,xt)D{^',x)di'{x) = j D{^,x)dv{x) j ,x)dr'{x) = ^ 

( 21 ) 

because linear combinations of the polynomials D{fi,x) are dense in Lfiivg). 
To prove (21) denote ^ T if the support of ^ and are disjoint, i.e., if 
there are no vertices i £ S which contain both particles from ^ and f,'. If 
^ T then under the measure vg, the polynomials D(^, •) and D{f^',■) are 
independent. Because of (20) it then follows, using duality and conservation 
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of the total number of particles in the dual process: 


lim f K^D{^,xt)D{^^,x)di^{x) 
t—^oo J 

= lim Vpt(^,C) / D{C,x)D{^',x)du{x) 

c 



= lim 

t —>-00 ^^ 

C-L?' 


= lim 

V:—^on ^ ^ 


= 0l«l+l?'l 


□ 


Notice that, for a single dual particle, that is to say when ^ = 5i, we 
have 

D{i,x) = ^. 

In the dual process, the motion of single dual particle is a continuous-time 
random walk jumping with rate p{i,j) from i to j. 

If we denote by Ptihj) the time t > 0 transition probability of this walk, 
then duality with a single dual particle implies the following “random walk” 
spread of the expected wealth at time t > 0. 


Proposition 3.2. In the model with generator (15), for all x ^ ^ and 
i G S we have 

E^{xi{t)) = '^pt{i,j)xj. 
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4 Generalized immediate exchange model 


Consider the update rule ([^ and assume that U and V are now independent 
and B{s,t) distributed (the original model is then recovered for s = t = 1). 

In other words, we consider the generator 

Ls,tf{x,y)= / / {fix{l-u) + yv,y{l-v) + xu)-f{x,y))(j)s,ti'^,v)dudv,{22) 


0 JO 


where 


<t>s,t{u,v) = 


B{s,t) 






(l-v) 


t-1 
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is the probability density of two independent B{s,t) distributed random 
variables. 

As before, the generator can be rewritten as L = P — I, where I is the 
identity and P the discrete Markov transition operator 

Ps,tfix,y) = „/ / / /(a:(l-n)+ ?/r;,y(l-u)+ xu)(/)^,t(u,u)d'udr;. 

B[s,ty Jo Jo 


In this generalized setting, the polynomials which we need for duality are 
now given by 


ds,t{k,x) 


x^T{s +1) 
r(s + t + A:) 


(23) 


and 


Ds,tin, m; x, y) = ds,t{n, x)ds,t{m, y). 


(24) 


These polynomials are associated to the Gamma distribution vp'^{dx) with 
shape parameter s + t, 


i,^+\dx) = 


r(s + t)0"+* 


dx 


via 

y ds,i(fc,x)r'g+*(dx) = 0*^. (25) 

As before, with a slight abuse of notation we also denote vp'^{dxdy) the 
product measure with marginals Vg'^^{dx). 

The same computation as the one following Q now yields that for a 
given k, I with 0 < k < n,0 < I < m, the dual process will jump from (n, m) 
towards {n — k + I, m — I + k), at rate 

( kl)— {k + s — l)l{m — I + t — l)l{n — k + t — l)l{l + s — 1)\ 

f's,t n,m, , B{s,ty {s + t + n — l)l{s + t + m — l)\{n — k)\kl{m — l)W. 

(26) 

where the factorials are to be interpreted as x! = r(x + 1), when x is non¬ 
integer. Notice that as before in Q we have that the rates sum up to one 

n m 

^ ^ m] k, 1) = 1. (27) 

fc =0 z=o 

This follows via rewriting 


with 


rs,t(n, m; k, 1) = Ws,t{n, k)ws,t{l, m) 


Ws,t{n,k) 


n\{k + s — l)\{n — k + t — 1)! 
B{s,t){s + t + n — V)^.k\{n — k)\ 
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and recognizing the probability mass function of the Beta binomial distri¬ 
bution with parameters n, s, t, given by 

BetaBin(re,s,t)(fe) = /(I - 

as a consequence one has 


'^Ws,tin,k) = 1 
k=0 


We can then state the generalized duality result, and its consequences, 


as in theorem 2.1 The dual process when initiated at (n,m) is once more 


an irreducible continuous-time Markov chain on the finite set {{k,l) : k + l = 
n + m} which converges to unique stationary distribution which we denote 
by denote v^^{k,l) and is given by 


<Xm{k,l) = 


r{s + t + k) T{s + t + 1) 1 

r(s + t)ki T{s + t)i\ 


where 


^s+t _ 
•^n+m 


E 


k^l‘.k-\-l=n-\-nn 


r(s + t + k) r(s + t + i) 
r(s-|-t)A;! r(s-|-t)/! 


(28) 


(29) 


Notice now that we have the analogue of (14), i.e., if we the product measure 
r'Q"’~*(fe, 1) conditioned on A; -|- / = n -|- m then 


D{k, l-,x, y)u^XLik, 1) = 


{x + y) 


n+m 


k,l\k+l=n-\-m 


{n + m)\^Xtln 


(30) 


is only a function oi x + y. As a consequence, we obtain the following result 
in the generalized model. 


Theorem 4.1. 1. The processes {Nt,Mt) and {Xt,Yt) with generator 

(22) and rates ( |26[ ) are dual with duality function (24). This means 
that, for all {n,m) G and {x,y) G [0, oo)^, we have 




S^t 

n,m 


Ds^t{Nt,Mt,x,y). 


2. As a consequence, the product measure ^^^{dxdy) is invariant. 

3. Moreover, starting from any initial state {x,y), the process (Xt,Yt) 
converges in distribution to {Se, ^(l — e)) where e is Beta{s + t,s + t)- 
distributed, and S = x + y = Xq + Yq. 

4 . The invariant measures with finite moments are of the form (Se, ^(l — 
e)), with e is Beta{s + t,s + t)-distributed. 
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We can then build the analogue of this model for many agents associated 
to the vertices of a graph S, as in equations (15) and (16). First notice that 
for a single dual particle, when ^ = 5i, we get 

D{i,x) = 


Xi 


s 1 


Just as before, the motion of single dual particle in the dual process is a 
continuous-time random walk, jumping with rate p{i,j) from i to j. If we 
denote by Ptihj) the time t > 0 transition probability of this walk, we then 
have the following result. 


Proposition 4.1. In the model with generator (15), for all x ^ It., i ^ S 
we have, for all r > 0 


K’\xi{r)) = Y^pr{i,j)xj{0). 
j 


5 Properties of the discrete dual process 

The discrete dual process is a redistribution model of independent interest. 
It is a discrete redistribution model of the same type as the original con¬ 
tinuous model in the spirit of the KMP process and its analogues of [T], 
where also a similar discrete process is the dual of the original continuous 
redistribution model. It is therefore useful to understand better the discrete 
dual process and its connection to the original process. 


5.1 Reversible measures 


Define the discrete Gamma distribution with shape parameter s+t and scale 
parameter 0 < 0 < 1 as the probability measure on N with probability mass 
function 




1 S" r(j + t + n) 
Zq nl r(s -|- t) 


(31) 


where Zq = (1 — 9) ^ ^ is the normalizing factor. We first recall that the 
dual process has generator 


f{n, m) = Ts,t{'n, m; k, 1) {f{n — k + l,m — I + k) — f{n, m)) 


k=01=0 


(32) 

where the rates are given by (26). It is important to notice here that this 
generator can be rewritten as follows 


.Sf f{n, m) = E/(n — Xi + X 2 , m — X 2 + Xi) — f{n, m) (33) 

(n) 

where Xi = Xj is Beta binomial distributed with parameters n, s, t and 
X 2 = x[^'^ independent Beta binomial with parameters m, s, t, and E is 
expectation w.r.t. these variables. 
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Proposition 5.1. For all 9 G (0,1), the product probability measures with 
marginals are reversible for the process with generator (32). 


Proof. The reversibility of Pg"*"* for the generator follows from a stan¬ 
dard detailed balance computation. Indeed, fix two configurations (n, m) 
and G with n -|- m = n' m!\ now, for any 0 < < n and 

0 < I < m such that n' = n — k + I and m' = m — I + k, it trivially follows 
that l<n' = n — k + l and k<m' = m — l + k and n = n' — I + k, 
m = m! — k + 1. In other words, for each redistribution of (n, m) according 
to r{n, m; k, 1), we can find a ’’reverse” redistribution of (n', m') according to 
r{n',m';l,k). Furthermore, these two redistributions are indeed reversible, 
as one may see by explicit computation, combining (p^ and (31) that 


r(n, m; k, l)i/n'^^{n)i'n'^^{m) = r{n+l—k, m+k—l; I, k)i'n'^^{n+l—k)i/n'^^{m+k—l) 


which implies detailed balance and thus reversibility. □ 


5.2 Scaling limit 

The fact that the rescaled Beta Binomial converges to the Beta distribution 
(by the law of large numbers) provides a connection between the discrete 
dual process and the continuous process. The continuous process arises as 
a limit of the discrete dual process where the number of initial “coins” is 
suitably rescaled to infinity. This is expressed in the following result. 

Theorem 5.1. Let nK,rnK be a sequence of integers indexed by K G N, 
and such that 

UK rriK 

as K ^ oo. Then we have that the corresponding processes nx(t) / K, mx{t)/K, 
with generator ( [3^ converge to the continuous process with generator ( p2) ), 
starting from (x, y). 

Proof. Define a number A > x + y. Because convergence of generators 
on a core implies convergence of the processes, it suffices to show that for 
smooth / : [0, ^4]^ —)• M 


lim {AL’fx){nx,rnx) = Ls,tf{x,y) (34) 

K^oo 


where fx{n, m) = f{n/K,m/K), is given by (32), and Lgg, by (22). Con¬ 
sider Beta binomial with parameters nx, s, t, and independent 

Beta binomial with parameters mx^ s, t. By the law of large numbers it fol¬ 
lows that 

^ xYs,t, -y;- ^ yYf t 


K 


K 
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with Yg^t, Ygt being independent distributed. Therefore, by smooth¬ 
ness of / and dominated convergence, as —)• oo we have 


lim E(/i^(nx - 
K^oo 


= E{f{x - xYg^t + yY'^t^y - yY^t + xYg^t) 
= Ls,tf{x,y) + f{x,y) 


which shows (34). 


□ 


6 Self duality and SU{1,1) symmetry of the dual 
process 


In this section we show self-duality with the self-duality polynomials which 
are naturally associated to the reversible discrete Gamma distributions. 
More precisely, we define the following discrete polynomials: 


ds,t{k,n) 


n\ r(s-|-t) 

(n — A:)! r(s + t + k) 


(35) 


where negative factorials are defined to be infinite. These polynomials are 
naturally connected to the discrete reversible Gamma distribution via 


^ dg^tik, n)z.^g+*(n) = p{0)^ (36) 


with p{9) = 6/{I — 9). Next we have the associated polynomial in two 
variables: 

Dg^t{k,l] n, m) = dg^t{k, n)ds^t{l, m) (37) 

Notice that in the case n = = [N'yJ, divided by and in the 

limit N ^ oo, these discrete polymials converge to the duality polynomials 
( |24[ ). We recall that the dual process has a generator of the form 

n m 

^f{n,m) = EE’'.. t{n,m-,k,l){f{n—k+l,m+k—l)—f{n,m)) = {P—I)f{n,m) 

k=01=0 

where the discrete transition operator 

n m 

Pf{n, m) = EE’'.. t{n, m; k, l)f{n — k + l,m + k — 1) 
k=0 1=0 

is indeed a Markov transition operator because, as we showed before, 

n m 

J2J2^^’t{n,m;k,l) = 1 . 

k=01=0 
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To prove self-duality of the process with generator (32), we show that it 
commutes with a SU{1, 1) raising operator + K 2 , from which we can 
generate the self-duality function via the strategy described in |2], namely 
by acting with on a cheap self-duality function coming from the 

reversible product measure. 

In order to proceed with this, we introduce the SU{1, 1) raising operators 

m, 

K'^f{n) = {s + t + n)f{n + l). (38) 

For a function f{n,m) of two discrete variables, we denote , resp. 
the operator defined in (38) working on the first (resp. second) variable. 
Similarly we have the lowering and diagonal operators 


K /(n) = n/(n - 1), f{n) = + n) f{n). 


(39) 


Together, the K , iF+, generate a discrete (left) representation of SU{1, 1); 
i.e. they satisfy the SU{1, 1) commutation relations 

[K+,K-] = 2K°, [K^, K°] = ±K^. (40) 

where [A, B] = AB — BA denotes the commutator. We will show in this 
subsection that the generator ^ defined in (32) has SU{1, 1) symmetry and 
that the self-duality follows as a consequence, in the spirit of [9], [I]. We 
start by noticing that by reversibility of the measure the function 


D{n', m!] n, m) = 5n',n5m'- 


n\T{s + t) m!r(s-|-t) 
''r(s + t + n) r(s + t + m) 


is a “cheap” self-duality function m, m Furthermore, we remark that the 
claimed self-duality polynomials can be obtained via 

D{n', m'] n, m) = "’"^2 T){n', m; n, m) 

where the operator ig working on the n', m! variables. Therefore, 

in order to prove that self-duality holds with the claimed polynomials (35), 
(37), it suffices to prove that commutes with the generator. In¬ 

deed, then from the general theory developed in [9], see also [2], it follows 
that e^i '*'^2 T){k,l',n,m), which arises from the action of a symmetry (an 
operator commuting with the generator) on a self-duality function, is again 
a self-duality function. 

Theorem 6.1. The generator in (32) and the operator Kf + K 2 com¬ 
mute, i.e., for a// / : —>• M we have 


.^{Kt + Kt)f = {Kt + 


(41) 
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Remark 6.1 (Hypergeometric Functions). We briefly recall some defini¬ 
tions and properties about hypergeometric functions we will need in the proof 


geometric function 2 F 1 
a b 


a b, 


of Theorem 6.1, On a suitable subdomain of {z € C : lR{z) > 0}, the hyper- 

’; z is defined via the following series expansion 

1 if k = 0 


2F1 


{a)k{b)k / N 

^ {c)k kV ^ 1 7 ^( 7 . _)_ 1 )... (^ 7 -_|__ 1 ) if k > 0 . 


k=0 


Note that for all n,k and t G 




{-n)k = (-l)^n • (n - 1 ) • • • (n - /c + 1 ) = (- 1 ) 


k r(n + 1) 
r(n — A: + 1) 


and 


(1 -n-t)k = (- 1 ) 


k r(n + t) 
r(n — k -\-t)' 


Moreover, as a particular case of Gauss’s summation theorem m, Theorem 
2 ), we can state that 


2F1 


—n s 




r(t)r(n + s + t) 
r(s + t)r(n + 1 ) ’ 


1 — n — t 

Some useful formulas are listed below: 

r(s + k) r{t k) _ 


n G N, s,t > 0. 


E 

k=0 


r(i + k) r(i + n — k) 

_ r(s)r(n + t) i-n)kis)k 1 r(s)r(n + t) ^^^ 


r(n + 1) 


k=0 


(1 - n - t)k k\ 


r(n + 1) 


—n s 
1 — n — F 


E 

A:=0 


r(/c + s) r(?T, — k + t) f 9i 


T{k + l)T{n-k + l) \e 2 


-k 


T{s)T{n + t) 

r(n + 1) 


2F1 


—n s O 2 
1 — n — t ’ 6*1 


s^T{k + s)T{n-k + t) fOA ^ , ^^^(n + ^-l) 


r(n) 

Proof. Let us prove that for all functions / : —>• M and (n, m) G 

P {K^ + K:^) f{n, m) = {K+ + K+) Pf{n, m). (42) 

By straightforward computations and substitutions, if we adopt the notation 

r(o + s + t) 


-n + 1 s + 1 _ 6*2 
2 — n — t ’ 6*1 


S.t 


T{b + s)r(a — 6 + t) ’ 


a > b > 0, s,t > 0, 


17 






























the l.h.s. rewrites {K'^ ;= 

n m 

PK^ f{n, m) = EE“V,(n.tK,,(m. 1 ) (+ K 2 ') {n — k + l,m — I + k) 
k=o 1=0 

^ n m 

TTi^EE 




k= 01=0 



n\ j m 

.k [l 


n 

L^J . L 
L J s.t L 




(s + t + (n — fc + l))f{n — k + I + l^m — I + k)+ 


+ {s + t + {m — I + k))f{n — k + l,m — I + k + 1) ], 


while the r.h.s 


K+Pf{n, m) = (s + t + n)Pf{n + 1, m) + (s + t + m)Pf{n, m + 1) 

n +1 m 


s + t + n I n + 1\ f m 


k' 



k’=0l'=0 

■ f{n + l-k' + + k') + 

n m+l 


n + 1 
k' 


S.t L 


S,t J 


+ 


s + t + m 


EE 


B(s,tP ^ ^ \k" 

'' ’ ^ k"=0l"=0 



n \ I m + 1 
I" 



n 

k” 


- S^t - 


m + l 
I" 


- 


■ f{n — k" + l", m + 1 — l” + k"), 

Let us introduce another shortcut: 

+ fcEN, s>0. 

As it is enough to show the identity only for the functions / : 
the form 


in 


f{n,m) ■.= 6 'i, 6»2 G (0,1), (n,m)EN^, 


we can recast (42) as follows: 

I I ^ 

n\m\ 


, I .ir/ , , .1 Yj Y - k)zs{l)zt{m - 1)- 

1 (^n + s + rjl (m + s + t) 

({s + t + {n-k + + {s + t + {m-i + k))e^-^+^e^-^+^e-^ = 


n +1 m 


r („^ E E -('=)-(-+1 - 


+ 


k=0 1=0 


k=0 1=0 
n m+l 
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EE Zs{k)zt{n - k)zs{l)zt{m - 1) 


k=01=0 



Oi 


{n + s+ t) - (n + l) 


02 

n — k + t 
n — k + 1 


l-k 


+ 02 


m + s + t — {m + 1 ) 


m — I + t 
m — I + 1 


+ 


+ {01 - 02) EE Zs{k)zt{n - k)zs{l)zt{m “ 0 ( 


l-k 


k=o /=o 


02 . 


0i{n + l)zs{n + l)2;t(0) ( ^ 


-(n+l) m 


Zs{l)zt{m - 1) ( ^ ) + 


1=0 
m+1 n 


02 . 


/a \ m^i /n 

+ 02{m + l)zs{m + l)zt{O) {^-^j '^Zs{k)zt{n-k) ^ 


k =0 



-k 


Since 


, ,, n — A: +1 , k 

n + s + t-{n + 1)- - -- = s - (l-t) 


n — k + 1 


n — /c + 1 ’ 


we can further simplify 

n m 


s{0l+02) EE Zs{k)zt{n - k)zs{l)zt{m - 1 ) 


k =0 1=0 
n m 



l-k 


+ 


(i-*>EE" s{k)zt{n - k)zs{l)zt{m - 1) 


k =0 1=0 

n m 


l-k 


+ {01 - 02) EE Zs{k)zt{n - k)zs{l)zt{m “ 0 ( 


01 


+ 02 - 


n — k + 1 m — I + 1 


l-k 


k= 01=0 


02 . 


0 i{n + l)zs{n + l) 2 ;t( 0 ) ( ^ 


-(n+l) m 


Yj Zs{l)zt{m - 1) ( ^ ) + 


1=0 
m+1 n 


02 . 


/0 \ /Q 

+ 02{m + l)zs{m + l)zt{O) '^Zs{k)zt{n-k) ^ 


Now, by noting that 
k 1 


and 


k =0 


1 



-k 


1 


r(A: + l) r(A:) r(n-A: + l)(n-A: + l) r(n - A: + 2) ’ 

and by using the shortcuts 

n f 0 \ ~^ ™ /o''’' 

N :='S^ Zs{k)zt{n — k) ( and M :='S^ Zs{l)zt{m — l) ^ 
k =0 ^ 1=0 
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- -^Tik + s) r{n - k+ t) f di\ ^ X ” r(/c + s) r(n - A:+ t) ^ 

■-E r(t) rin-t + ijis^j “ ^-E rw r{n-j: + 2) ie.J 


and similarly for M and M, we can continue with 

r(n + i + s)r(t )/01 


M sOiN + (1 - t)eiN - (01 - 02)A^ - 02- 


r(n + 1) 


9, 


(43) 


= N -S02M - (1 - t) 62 M - (01 - 02)M + 01 


r(m +1 + s)r(t) /01 


r(m +1) 


02 


Note that, as in Remark 6.1, we can rewrite these quantities N, N etc., in 
terms of hypergeometric functions as follows 

—n s 02 

^(71 + 1)^”^^ 1 —n —t’0i 






1 — n 1 + s 02 
2 — n — t ’01 


and 


N = 


02 1 


r(s + i)r(n -1 + t) 2 Fi 


0ir(n + l) 

Therefore, the expression 

s0iiV + (1 - t)eiN - (01 - 02)iV - 02 


—n 1 + s 02 
2 — n — f ’ 01 


- r(t - l)r(n + 1 + s) 


r(n + i + s)r(t) /01 


r(n + l) V02 


simplifies to 

r(n + f — 1) 


r(s + i)- 


r(n + 1) 

+ 02(1 — t)2Fi 


9i{n + t- 1 ) 2^1 

—n 1 + s 02 
2 — n — t' 9i 


—n s 02 
1 — n — t ’ 01 
02 


+ 


- 02n ( 1 - — ) 2 T 1 


1 — n 1 + s 02 
2 — n — t ’01 
(44) 


—n 1 + s 02 

-n fl - 2 T 1 

1 — n 1 + s 02 

1 — n — t ’ 01 

V 0 iy 

2 — n — t ’01 


By some standard manipulations of hypergeometric functions, the expression 
(1 — t) 2 Fi 
reduces to 

/ , ^ r —n s 02 

-{n + t-lhF, 
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In conclusion, if we go back and plug the latter expression into (44), we can 
rewrite the l.h.s. in (43) as 


r{m + t) 


r(m + 1 ) ■ 


—n s 02 
1 — n — t’ 01 


r(s + i) 


r(n + t — 1) 
r(n + 1) 


{01 - 02 )in + t- 1 ) 2 Fi 


—n s 02 
1 — n — t' 01 


{01 - 02)sT{s) 


2 r(n + t) V{m + t) 
r(n + 1) r(m + 1 )' 


:Fi 


—n s 02 

2 F 1 

—m s 

0 _i 

1 — n — t' 01 

1 — m — t 

02 _ 


By simply replacing n by m, 0i by 02 etc. and exchanging the sign in the 
latter expression, one simply obtains the explicit form of the r.h.s. in (43), 
which indeed proves identity (42). □ 


We extend now the commntation of the generator with + K 2 to full 
SU{1, 1) symmetry of both the discrete and the continuous model. For this 
we need some additional notation. Denoting the operators working on 
functions / : [ 0 , 00 ) —>■ M 


jr+/(x) 

= xf{x) 

(45) 

f{x) 

= {xdl + {s + t)d,r;)f{x) 

(46) 

JF^f{x) 

= (x+^)/(x) 

(47) 


we have that the algebra generated by forms a (right) representation 


of SU{1, 1), i.e., satisfy the commutation relations (40) with opposite sign. 
Moreover, this continnons right representation is linked with the discrete 


left representation used before via the duality polynomials (23), i.e. 


JF°'ds,t{n, x) = K°‘ds,t{n, x), a G {+, -, 0 } 


(48) 


where works on x, and K on n (see e.g. [2] for the proof). 

We now first formulate a simple lemma, showing that 0~^K~ is the 
adjoint of in L‘^{vb). 


Lemma 6.1. 


Let 


he the reversible measure for the discrete dual process, 


defined in (31). We have in L'^{iXq 


s+£\ 


{K*r = -^K- 


where K°‘ are the operators introduced in (38), (39) 


Proof. Let /, <7 : N —M be functions with compact support, then we 
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compute 


'^f{n)K+g{n)iy"g'^\n) 

n>0 


1 


f{n){n + s + t)g{n + 

n>0 


9^r{s + t + n) 
n! r{s + t) 


1 


'^f{n)g{n + l) 

n>0 


9 '^ r(s + t + n + 1) 
n! r(s + t) 


1 1 

9Ye 


^re/(n - 

n>l 


l)g{n) 


9'^T{s + t + n) 
n\ T{s + t) 


= f{n)g{n)u^g-^\n) 

^ n>0 


□ 


We are now ready to prove the full SU{1, 1) symmetry of both the orig¬ 
inal continuous process and the discrete dual process. To explain this, we 
denote the coproduct 


A : ^{SU{1, 1)) ^ ^{SU{1, 1)) ® ^{SU{1, 1)) 


which is defined on the generators as A{K°^) := Kf + and extended 
to the algebra as a homomorphism. We then say that the process with 
generator L has full SU{1, 1) symmetry if it commutes with every element 
of the form A(A), A G ^{SU{1, 1)). This in turn follows if it holds for the 
generators iiT", by the bilinearity of the commutator. 


Theorem 6.2. Let denote the generator of the discrete dual process, 
defined in (32), and L the generator of the continuous process defined in 
(22). Then we have for a G {+, —,0} the commutation properties 

Kf + Kf] = [L, jxrfi + jefi] = 0 (49) 


As a consequence both ^ and L have full SU{1, 1) symmetry. 


Proof. We start with the discrete process. Because the sum of the wealths 
is conserved, Af trivially commutes with + K^- We showed in (41) that 
it commutes with ■ To show that it commutes with Kf + Kf 

we use lemma 6.1 and the fact that Af is self-adjoint in by the 

reversibility of Vg- 


[.if, Kf + Kf] = 9[^\ {Kt + K+T] = -9 (lAf, (K+ + K+)]y = 0 


We then turn to the continuous model, using ( [4^ . We show the commuta¬ 
tion with Jiry , the other cases are similar. We consider Dg^n, m; x, y), 
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the duality polynomial defined in (24), (25), and abbreviate it simply by D, 
where in what follows we tacitly understand that operators of the form 
are working on x, y and of the form K on n, m. In this notation, remark that 
operators working on different variables always commute (e.g. commutes 
with etc.). We can then proceed as follows, using duality which reads 
= LD. 




+ je'+)LD 


On the other hand, via ( |48[ ) 


L{K+ + K+)D 
{Kt+K+)^D 
^{Kt + K+)D 
^(jr+ + je'2^)D 


where in the third equality we used the commutation of ^ with . 

Combination of these computations then gives indeed 




on the functions D, and then by standard arguments on all / in 

□ 
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